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Abstract 

o ' 

The dilaton-limit fixed point and the scaling properties of hadrons in the close vicin- 
ity of the fixed point in dense baryonic matter uncovered in hidden local symmetry im- 
plemented with spontaneously broken scale symmetry are shown to reveal a surprisingly 
intricate interplay, hitherto unsuspected, between the origin of the bulk of proton mass and 
the renormalization- group flow of the w-nuclear interactions. This rends a theoretical sup- 
port to the previous (phenomenologically) observed correlation between the dropping nucleon 
| mass and the behavior of the w-nuclear interactions in dense matter described in terms of 

■ half skyrmions that appear at a density denoted 77,1/2 i n skyrmion crystals. The role of the 

, L(j-meson degree of freedom in the source for nucleon mass observed in this paper is highly 

reminiscent of its important role in the skyrmion description of nucleon mass in hidden local 
symmetric theory. One of the most notable novel results found in this paper is that the 
nucleon mass in dense baryonic medium undergoes a drop roughly linear in density up to a 
^ \ density (denoted h) slightly above nuclear matter density (no) and then stays more or less 

constant up to the dilaton limit fixed point. The possibility that we entertain is that n coin- 
cides with or at least close to 77,1/2. We note that this feature can be economically captured by 
the parity-doublet model for nucleons with the chiral-invariant mass mo ~ (0.7 — 0.8)m7v- It 
is found in one-loop renormalization-group analysis with the Lagrangian adopted that while 
the p-NN coupling "runs" in density, the cj-NN coupling does not scale: it will scale at two- 
loop or higher-loop order, but at a slower pace, so it is more appropriate to say it "walks" 
rather than runs. The former implies a drastic change in the nuclear tensor forces, affecting, 
among others, the nuclear symmetry energy and the latter generating the stiffness of the EoS 
at density higher than that of normal nuclear matter. 
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1 Introduction and Conclusion 



If hidden local symmetry (HLS for short) is assumed to hold in the vicinity of chiral restoration 
in dense baryonic matter and if U(2) symmetry is a good flavor symmetry for the vector 
mesons p and uj in medium as it is in free space, it follows from RG (renormalization-group) 
considerations that the U(2) hidden gauge coupling constant g and the "effective" V (vector- 
meson) -nucleon coupling constant qvnn will scale in density as the quark condensate (qq) 
scales, and consequently as (qq) — > in the chiral limit, both the p-nucleon and w-nucleon 
couplings as well as their masses will go to zero pp. This implies two dramatic effects in 
nuclear processes even slightly above nuclear matter density: (1) The suppression of the p- 
nuclear coupling will remove the p -exchange tensor force and hence strongly affect the nuclear 
symmetry energy [2]; (2) the suppression of the w-nucleon coupling will soften short-range 
repulsion in nuclear interactions and hence make the baryonic matter collapse at moderate 
density [3J. The first effect, when treated appropriately, turns out to be consistent with the 
EoS of compact-star matter that involves central densities of ~ (5 — 6)no (where tiq is the 
normal nuclear matter density), in fact playing a crucial role for explaining, within the HLS 
framework, the recently discovered 2-solar mass neutron star [4]. On the other hand, the 
second effect, if unmodified, would be disastrous for the stability of baryonic matter in the 
density regime relevant to compact stars. 

It is the purpose of this paper to suggest how to avoid the disaster due to the second 
(negative) effect without affecting the first (positive) effect. The key observation is that the 
U(2) symmetry, seemingly good in the matter-free vacuum, must break down in medium, and 
hence the properties of the isovector and isoscalar vector mesons behave markedly differently 
as density is increased. The key element in this is the origin of the bulk of proton mass 
that appears to have no direct link to chiral symmetry, its RG flow and its unsuspected 
association with the property of the w-NN interaction in dense medium. Our conclusion is 
that while the p and uj masses tend to zero (in the chiral limit), perhaps not in the same way, 
as density is increased, the effective w-nuclear coupling "walks" in contrast to the effective 
p-nuclear coupling that "runs." This feature was indicated in the phenomenology of compact- 
star matter studied in [3], and we show in this paper how that feature can be understood in 
the framework of HLS, e.g., vector manifestation and dilaton-limit fixed point. 

The basic assumption that we make is that local field theory can be applied to dense 
baryonic matter up to the density relevant to the EoS of compact stars. We will not, however, 
require that it be valid all the way to the chiral transition density denoted n c . In fact, we 
will not address what happens precisely at n c but consider approaching it from below. This 
means that we will not be able to properly account for the possibility of explicit quark degrees 
of freedom in discussing the EoS. 

We consider HLS Lagrangian that contains as relevant degrees of freedom, the pions and 
the lowest-lying vector mesons, p and uj. It may very well be that for realistic treatment, as 
stressed recently [5] , the infinite tower of hidden local symmetric vector mesons as indicated 
in holographic QCD models [6] need to be incorporated and the nucleons should be generated 
from such a generalized HLS Lagrangian. In this work, we will take the simplified Lagrangian 
in which the infinite-tower is integrated out leaving only the lowest vector mesons V = (p,uj) 
- in addion to pions, add baryon fields coupled to the mesons a la HLS and implement a 
scalar dilaton field x reflecting spontaneously broken scale symmetry linked to the QCD 
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trace anomaly. To be precise in definition, we shall call the dilaton-implemented HLS 
Lagrangian with baryon fields "dBHLS" while "(£BHLS" will stand for dilatonless baryon 
HLS Lagrangian. "HLS" will stand for the generic notion of hidden local symmetry as well 
as for meson-only theory (without dilaton). 

Our objective is to access nuclear matter as well as denser baryonic matter with a single 
Lagrangian, i.e., dBHLS. At zero density, that is, in the matter-free space, ^BHLS is "gauge 
equivalent" to baryon nonlinear sigma model and can be formulated, with the inclusion of 
chiral loops, to give a consistent chiral perturbation theory with baryons and vector mesons 
in a way parallel to HLS pQ. For this, scalar excitations, in principle, can be generated by 
loops. The QCD scalar qq is high-lying and hence does not figure in nonlinear sigma model 
in hadron dynamics at low-energy. Scalar glueball excitations will figure for the QCD trace 
anomaly, but they are also too massive. Thus the role of the dialton x m the vacuum is 
unclear and remains an unsolved problem. At tree order, however, the ^BHLS encodes the 
current algebra, and the dilaton x could be suitably interpreted in dBHLS as the lowest scalar 
excitation in low-energy pion dynamics. 

In going to nuclear matter density, one possible approach could be to do in-medium 
chiral perturbation theory with (£BHLS as one does with non-linear sigma model. However 
this requires high-order loop calculations, and this has not been done yet. What one can do 
instead is to do the mean field calculation with dBHLS as suggested in [8l [9] . The rationale 
there is that doing the mean-field with a chiral Lagrangian of the dBHLS type near the nuclear 
saturation density is equivalent to doing Landau Fermi-liquid fixed point theory provided the 
parameters of the Lagrangian are suitably scaled [HI EJ [TO] . What is required is that near the 
Fermi-liquid fixed point, the scalar x should be (predominantly) a chiral scalar with a mass 
around 600 MeV. 

Now given the dBHLS Lagrangian with the parameters sliding with density, normalized 
at the nuclear saturation density no, the question then is how to go to higher density close to 
the chiral transition point n c ? This question was raised and answered in [3]. It involves what 
is called "dilaton limit fixed point" (DLFP for short) - defined precisely later - introduced 
by Beane and van Kolck (TTJ. The basic idea is as follows: While the dilaton x should be 
dominated by chiral singlet component at low density near no in a complicated configuration 
consisting of multiquark and glueballs, as one approaches chiral restoration at n c , one should 
recover a Gell-Mann-Levy (GML)-type SU(2)l x 5t/(2)^linear sigma model with a scalar a 
making up the fourth component of the chiral four- vector (iri, ^2,^3, cr) |12j . 

For what follows, we need to generalize the GML sigma model to a parity-doubled model 
with nucleons in mirror assignment. In the original GML model (with nucleons in standard 
assignment), spontaneous breaking of chiral symmetry (in the chiral limit) entirely generates 
both the mass of the scalar a and the mass of the nucleon. Thus in the limit that the 
symmetry is restored, both the a and the nucleon should become massless. We will argue 
that while the a mass could go to zero to join the triplet pions, the nucleon mass need not 
be zero at the symmetry restoration. This feature can be captured by introducing parity- 
doubling in the baryon configurations. We shall call the former "standard" baryons and 

The role of scalar fields in effective Lagrangians is highly problematic in general and it is not at all obvious 
how to do this also in our case. We will however be guided by phenomenology in low-energy nuclear physics, 
namely, the EFT Lagrangian be treated at mean field with the parameters of the Lagrangian "sliding" with 
the density of the background in the spirit defined in [Tj- This will be the guiding principle in what follows. 
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the latter "parity-doubled" baryons. As will be elaborated in detail below, within the EFT 
framework, the nucleon mass cannot decrease much as density increases without getting into 
conflict with nature. To be specific, we write the nucleon mass parameter in the Lagrangian 

as 

7727V = mo + Am(K) (1) 

where m,Q is a chirally invariant mass term, a constant independent of the chiral condensate 
k oc (qq), and A is dynamically generated mass that tends to zero as k — > 0. The mass 
term uiq can appear in chiral Lagrangian without upsetting chiral symmetry provided parity- 
doublets are introduced [13]. A substantial tuq, indicated in nuclear phenomenology [I], will 
be the key issue in this paper. 

The dilaton limit fixed point (DLFP) arrived at in (H [14] in the parity-doublet baryon 
model dBHLS can be summarized as follows. As shown in [3] - and will be recalled in detail 
below, the idea is to make a field re-parametrization so that dBHLS Lagrangian with chiral 
symmetry in nonlinear realization can be linearized in the limit k — > 0. If one assumes that 
U(2) symmetry holds for p and u>, one arrives at the DLFP of the form 

(gv — 1, gA — L Am = m/v — mo) — > (0, 0, 0) as k — > (2) 

where V = (p,co), gA is the axial vector coupling constant and gv is the "induced" 1^-nucleon 
coupling defined by the effective vector-meson-nucleon coupling 

9vnn = g(gv - i) (3) 

where g is the U(2) hidden gauge coupling constant. As k — > 0, we expect that, even for 

gpNN = g^NN 0. (4) 

The "vector manifestation" (VM) of HLS (and (IBHLS) corresponds to g — > 0. We can see 
from ([3]) and (j4]) that the pNN coupling g p NN can go to zero before the VM/HLS fixed point 
(identified with chiral restoration) is reached. 

Several important consequences follow from the property (H|). 

If dBHLS were applied in the mean-field approximation to baryonic matter much denser 
than nuclear matter, then one would expect that the nuclear symmetry energy that figures 
importantly in the EoS for compact-star matter, proportional to g 2 p NN-> decrease and vanishe 
at the DLFP. One should however note that there is no solid argument why the mean-field 
approximation should hold at density much higher than the saturation density, i.e., Fermi- 
liquid fixed point, where mean-field approximation is valid. In fact, g p NN would imply 
that the p contribution to the nuclear tensor forces get suppressed, leaving the pion tensor 
force more effective. This feature turns out to lead to the stiffening - instead of the softening 
- of the symmetry energy, contrary to the mean-field expectation, as has been observed in 
effective field theory approach to nuclear dynamics. Indeed what was found in [3] provides a 
rather strong support for the prediction (j4j) for the p meson. These matters are discussed in 
SecEH 



4 



The situation with the w-NN coupling is quite different. As noted in [3], would 
imply the suppression of the principal mechanism in EFT involving oj exchanges for the 
short distance repulsion in nuclear interactions indispensable at high density. Since one 
cannot pinpoint the density at which the DLFP sets in, one cannot say that the predicted 
suppression of the w-exchange repulsion is inconsistent with nature. However it is generally 
believed in nuclear community that any significant reduction of the cj-NN coupling would 
make the EoS of nuclear matter too soft, causing difficulty in getting correct saturation. 
In [3], it was observed that unless the nucleon mass is suitably reduced, thereby increasing 
repulsion, it would be difficult to reduce the w-NN coupling. This suggests a close correlation 
between the behavior of the nucleon mass and the w-NN coupling in dense medium. This 
issue will be further elaborated on in Sec. 15.21 One could investigate the interplay between 
the two phenomenologically in the EFT formalism employed in [3]. In this paper, we will 
present a theoretical reasoning as to how the interplay can take place. 

Taking the hint that the U{2) symmetry can be substantially broken in medium, what 
we propose is to depart from hidden local symmetry in \U{2)i x i7(2)ft] g i bai X (2)v]iocai 
and consider [SU(2)l x 5^7(2)^] g i b a i x [SU{2)y x f7(l)v]i OC al- Denoting gauge couplings 
by g p and g u and induced couplings by gy p and gvu, an analysis parallel to the [U(2)]i oca i 
case [3j [14] reveals that as k — > 0, one approaches the DLFP 

(gv P ~ 1,9A - 9Vp,Am) -> (0,0,0), (5) 

which is what was found in [SHU]- But there is a major difference for the gy u — 1 coupling: 
There is no constraint for that coupling in going toward the DLFP. Furthermore as will be 
shown below, one-loop RG analysis shows that gy u does not "run" in contrast to gy p which 
drops rapidly to zero toward the DLFP. This leads to our conclusion that as density increases 
above no, the U{2) symmetry must be broken down significantly, and the g^NN coupling does 
not drop as fast as g p NN does. 

In short, the RG analysis with ^BHLS and the mean field treatment with dBHLS, both 
given in this paper, the large N c analysis of dense skyrmion matter described in |15| and 
the phenomenological study of the EoS for compact-star matter of [1] all converge to the 
conclusion that the dropping of the in- medium nucleon mass stops at h « ni/ 2 ~ 2riQ with 
the mass staying constant up to near chiral restoration, suggesting an itiq ~ (0.7— 0.8)m7v in 
the parity-doublet model for baryons. This property is found to have an intimate connection 
with the role that the cj-meson degree of freedom plays in the structure of nucleon and nuclear 
matter. 

In what follows, we provide details to what are given in sketch above. 

2 Hidden Local Symmetric Parity-Doublet Model 

In this section, we give a precise definition of the model we will study. We will first write 
down and discuss the hidden local symmetric parity doublet model without dilaton field, i.e., 
(iBHLS |14j . We will be focusing on the case where the chiral invariant mass tjiq in (|26p is 
non- vanishing, in fact, substantially big, but we will also discuss the case for uiq = 0, i.e., 
"standard" scenario. 
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We will limit our considerations to two flavors (Nf = 2). Motivated by the finding in [3], 
we relax the U(2)y symmetry for p and u mesons. The symmetry of the Lagangian involved 
is G g iobai x #iocai, where G g i oba i = [SU(2) L x SU (2) R] g i oha \ is the global chiral symmetry and 
-ffiocal = [SU(2)v x U(l)v]\oca,\ is the hidden local symmetry. We take p as the gauge bosons 
of [SU (2) v]\ocai and w as of [£7(l)v]iocai- The basic quantities are the HLS gauge bosons, Vp 
and Vjt, 



9 P P' 



T 

2 



(6) 
(7) 



where g p and are the HLS gauge couplings that will be taken unequal for the local SU(2) x 
U(l) symmetry concerned. The vectors transform 



h(Tr(x),g L ,9R) G [SU(2) v ] local and u = u(ir(x), g L , g R ) G [^(l)v]i oc 



(8) 
(9) 



with h = h(Tr(x),g L ,g R ) G [SU(2) v ] local and u = u(tt(x), g L , g B ) G [U(l)v]\ oca x in terms of 
9l,r G [*S'^ r (2)-L,-R]gi bai an< ^ the two matrix valued variables £l and £r, combined in a 2 x 2 
special-unitary matrix representing the pion field 



U 



transforming 

The variables £'s transform as 

They may be parameterized as 

Cl,r 



t 



(10) 

(11) 

(12) 
(13) 



Here 7r = vr a; y denote the pseudoscalar Nambu-Goldstone (NG) bosons associated with the 
spontaneous breaking of G g i b a i chiral symmetry, and a^ and a p (= ap 11 ^) are the NG bosons 
associated with the spontaneous symmetry breaking of J7(l)i OC ai and <S?7(2)i oca i respectively. 
The as are absorbed into the HLS gauge bosons through the Higgs mechanism, giving rise 
to their HLS boson masses. F n , F atJ and F up are decay constants of the associated particles. 

To construct hidden local symmetric Lagrangian, it is convenient to introduce the Maurer- 
C art an 1-forms 



a 



~ p 



1 

2i 
l r 



2/ 



d^r ■ 4 + d»zl ■ d 



with the covariant derivatives of £l r 



D^r = d^ R - %Vtt R - iVfcn . 



(14) 



(15) 
(16) 
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transforming homogeneously, 

d^ .I — > uh d^ || W v) , (17) 

and a± £ G — H and dn G 

With the above definitions, we can immediately write down the Lagrangian for the mesonic 
sector. To the leading order in derivative expansion (i.e., to 0(p 2 )), it is 



C M = F%tv [d_L M d^] + F 2 p tr a\\^ + F 2 tr [dyj tr 



^tr [V^l " T2 iT tr t^H . ( 18 ) 



with 



7?2 _ p2 







2 

1 / 1 



So 2 Vfi£ 9p 



(19) 



and the field strength is given by 

ynu = d pyv _ fftyn _ i [yp 5 , (20) 
y = + y£ . (21) 

We recover local U{2) symmetry if we set F 2 U = F 2 p and g u = g p = g in (fTKj) . 

One can readily write down, following [31 [14], the Lagrangian for parity-doublet nucleons 
coupled to HLS vectors 

C N = Qvy^DfjQ - giF n QQ + g 2 F n Qp 3 Q - im Qp2j 5 Q 

+9VpQ^a llfl Q + gvoQ^ti [a^] Q + g A Qp^a ± ^ b Q , (22) 

with the nucleon doublet Q = ( J transforming 



Q2 

Q^uhQ, (23) 

where the covariant derivative of Q is 

D»Q = (dp - iV P p - iV ull ) Q, (24) 

and the pi are the Pauli matrices acting on the parity-doublet. gA, gvp and gvo = \(gvu — gvp) 
are dimensionless parameters. To diagonalize the mass term in Eq. (|22p . we transform Q into 
a new field N: 

* + \- 1 ( ^ ^~^\( Ql \ 

N- ) ~ ^sh5 V T,e^ 2 -e 5 / 2 ){ Q 2 J ' (25) 

where sinh<5 = We identify N± as parity-even and parity-odd states respectively. The 

nucleon masses are found to be 



m N± = Tft^ + y{giF 7r y + m 2 , (26) 

co S h6= mN + +mN - . (27) 
2m V ' 
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Finally, we arrive at the Lagrangian in parity eigenstate as 

C N = NipN - NMN + g A N^Ga ±fM -f 5 N 
+ g Vp N^a^N + gvoNjHr [ay] N , 
m>N + \ A / tanh 5 75 / cosh 5 



M 



rriN_ 



G 



75 / cosh S — tanh 5 



(28) 
(29) 



It is convenient for the analysis of u>-nucleon coupling to change slightly the Lagrangian 
28l) . We define quantities belonging to the algebra of G g i b a i X [SU (2)v] local a s 



off = 1 - 
" 2i 



where 



(30) 
(31) 

(32) 



and 



a and a are related to each other via 



(33) 



a 1 



o 1 



6"' 



<7u 



2F n 



Then, the Lagrangians (|18p and (|28[) take the form 



-Al 



F> [a ± ^ ± ] + F CT 2 p tr 



Q!|UQ!|7 



+ I -£ 3<^ M 



(34) 
(35) 



(36) 



-AT 



iVi^iV - NMN + g A N^Ga ±ll ~f 5 N 



+ 9VpN'fa\\ vt N + g Vu Ni l 
where gvo is replaced by ^ (gvw — gvp) and 



^T 1 



d»p v ■ \ - d v p> x ■ \ 



U 2 



r 2 ' " 2 



(37) 



(38) 
(39) 



Note that the u meson couples to nucleon and NG cr^, but there is no u coupling to other 
mesons, i.e., tt and p, at tree order. There can be tree-order W7r 3 and uj-ir-p couplings in the 
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homogeneous Wess-Zumino term in the anomalous part of the HLS Lagrangian that could 
give rise to a one-loop correction to w-nucleon coupling, but does not contribute at the order 
we are working with. We willl see in Seed] that at the one-loop order, the coupling gy u does 
not scale. 

One can read off the vector meson mass and the pmr coupling constant at tree level as 

m 2 p = g 2 p F 2 p , ml=glF^, (40) 
gp-K-n = 7, a p9p > ( 41 ) 



F 2 
— 12. 

a P ~ 7?2 

7T 

The vector mesons couple to nucleons as 



(42) 



9 p n + n + = g P N-N- = (gvp — i)g P , (43) 

gu)N + N + = gu)N_N_ = {gvu — i) g^i (44) 

and the axial vectors coupling as 

gAN+N+ = -3ajv_jv_ = ffAtanh<5 (45) 

where the subscripts stand for the parity of the nucleon doublet. When U(2) symmetry 
is restored, we will have g^NN = gpNN and the w-nuclear coupling will vanish as in [3] in 
approaching the DLFP. In what follows, we will not assume U(2) symmetry. 



3 Going Towards the Dilaton-Limit Fixed Point 



In order to study what happens to the baryonic matter as density increases, we need to incor- 
porate the dilaton field x that represents spontaneously broken scale symmetry of QCD. The 
explicit scale symmetry breaking associated with the trace anomaly that is also presumably 
responsible for the spontaneous breaking [15] will not figure directly in our consideration. 

We follow the standard trick of inserting the "conformal compensator" field x i n t° the 
Lagrangian (I18D and (I22p to obtain scale symmetric Lagrangian, with the scale invariance 
broken spontaneously. It is given by 



C 



C 



M 



-X ' 



X 



X 



Qij^DnQ - giF n —QQ + giF^—QpzQ - im Qp 2 ^Q 



r x r x 



+gv P Qi ll 6i\\pQ + gvoQi^r [&y] Q + gAQp3i^a±^ 5 Q , 



F 2 F 2 
-^X 2 tr [a^&l] +-^-X 2 tr 

x 



\ 



a \\v \\ 



+ 2F 2 X tr [, 



an^J tr 



Q i 



1 



1 



tr [uj^ v 



;dpX • d^ X ~ V(x) 



(46) 
(47) 

(48) 
(49) 
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where V(x) is the Coleman- Weinberg- type dilaton potential that breaks scale symmetry 
spontaneously. We do not write down its explicit form since it is not needed for our purpose. 
Here, F x is the vacuum expectation value of x a t zero temperature and density. 

To move towards a chiral symmetric GML-type linear sigma model, we do the field re- 
parametrizations E = f/xp 5 " = s + if ■ tt - that also defines scalar s - and 



or equivalently 



W 

Q 



1,2 - 7T 



d+ei)±7 5 (4-£i)]Qi,2: 



d + A) + ^375 (4 - el 



2 

2 [(& + & + p3i5 (£ R - a)] ^ • 



(50) 

(51) 
(52) 



With these reparametrized fields and going to parity eigenstates, one finds a complicated 
expression for (|4"6|) composed of a part that is regular, £ rcg , and a part that is singular, £ s ing, 
as tr(EE') = k 2 = 2 (s 2 + 7r a2 ) — >■ 0, where a is iso-spin index. The singular part that arises 
solely from the scale invariant part of the original Lagrangian ()46|) has the form 



Aing = (gv P ~ 9a)A[ 1/tr 



EE 1 " 



+ (a - 1)8 ( 1/tr 



EE 1 " 



(53) 



p2 

where a = and 

X 

A 



to 



tr (x^ 1 ")) {e, E 1 "} - 2tr (EE 1 ") (s^E 1 " + E^E 



-tr 



^375 (s^St - E^Ej ^ 
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1 



16a 



tr ( EE f ] tr 



(EE* 



tr 



[ EE f 



(54) 
(55) 



That /^sing be absent leads to the conditions that 



9V P ~ 9 A -> , a - 1 -)• . (56) 

Using large N c sum-rule arguments |11| and the RGE given in the next section , we infeJl 

g A -l^0. (57) 

In the density regime where GML-type linear sigma model is valid, the nucleon mass can be 
given as 

(58) 



mN ± = T92(s) + y (gi{s)) + m 



'0 • 



2 It is most plausible that qa — 1 = is reached only after qa — gv P — is reached. This is because 
tr(EE') — y gives the constraint qa — gv P = 0, but not gA — 1 = 0. However the RGE given in the next 
section has the infrared fixed point gv P — 1 = 0, so the point at which gA — gv P = could be very near the 
chiral restoration point at which gA is close to 1. In fact it is observed phenomenologically in Gamow- Teller 
transitions in heavy nuclei that jJsjI. 
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where (s) is the vacuum expectation value of s. As the chiral symmetry restoration point is 
approached, (s) — > 0, so in the limit tr(EI^) — > 0, we expect 

vn 7y_j_ — ^ mo . (59) 

These are the constraints that lead to the dilaton limit as in [3] and announced above. It 
follows then that 

g P NN = g P (gv P - l) ->■ 0. (60) 

We thus find that in the dilaton limit, the p meson decouples from the nucleon. In contrast, 
the limiting tr(SS^) — > does not give any constraint on (gvuj — 1)- The w-nucleon coupling 
remains non-vanishing in the Lagrangian which in unitary gauge (with cr u = a p = 0) and in 
terms of fluctuations s and tt around their expectation values, takes the form 

C N = Ni@N - NMN - gi N (ds + p 3 j 5 iT ■ Pj N 

+ 92M ( P3 s + G l5 ir ■ N + (l - g V u) gujN~^N , (61) 

which is the same as the Lagrangian given in [13] except for the w-nucleon interaction. This 
is just the nucleon part of the linear sigma model in which the ui is minimally coupled to the 
nucleon, applicable infinitesimally below the critical density n c with the mass M replacing 
mo- What is significant with this result is that it shows that the suppression of the w-repulsion 
predicted with £7(2) could be absent due to potentially significant U (2) symmetry breaking 
as is indicated in the phenomenology of EoS for dense baryonic matter. We will return to 
these matters in Sees. 15.2 1 and 15.31 



4 RG Analysis of the cj-Nucleon Coupling 

We have arrived at the DLFP ([2]) by linearizing dBHLS Lagrangian (j46p (which is gauge 
equivalent to nonlinear sigma model) in the mean-field approximation. Very near the fixed 
point, the resulting effective Lagrangian is a GML-type linear sigma model, to which the uj 
meson is minimally coupled. The isovector vector meson p is decoupled in the limit. 

In this section, we interpret the DLFP in terms of RG (renormalization-group) flow. To 
do this, we take the hidden local symmetry Lagrangian with baryons but without the dilaton 
X (i.e., (ffiHLS). As mentioned, in the chiral perturbation approach that we will be taking 
in this section, the role of scalar dilaton is problematic if introduced naively. We will first 
discuss the standard assignment for the nucleon, which gives a clearer picture of what's going 
on, and then consider the mirror assignment with parity-doubling. 

4.1 RGEs in the Standard Assignment 

In the standard assignment with rriQ = 0, the nucleon and its parity partner are decoupled. We 
shall consider only the positive-parity nucleon (set ./V = N + ) in analyzing the RG properties 
of its coupling to mesons. The calculation is straightforward, the only difference from what 
was done in [TJ] being that we have SU(2) x U(l) local symmetry instead of U(2). Omitting 
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the details that are given in |14| - apart from terms involving the uj meson, we simply write 
down the RGEs at one-loop order, 



drriN SrriN p 

^ dfl 167T 2 " 



(62) 



dgA 



m 



N 



8vr 2 F 2 



•Fo + (l 



9Vp) ^2 Fi + (g Vp - 9a) i~2 F 2 , 



8tt 2 " 



8vr 2 " 

9 A n 



d9A m" N g A ~ g A ~ . 2 , g A ~ 



(63) 
(64) 



Mt- {gvuj - i) 
djjL 



(65) 



The explicit expressions of JFi and Qi are given by 



•F m 
F 
Fl 

Gi 
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^ (u 2 
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1 
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+ 



2^ 2 P 



3 K,^ 

2 I ^ 2 ^ 2 P 



m JV + ffp ( 2 - 



*" ^ 2 



2F 2 



1 



9v P 



->g A + --^ + 2gv P 



Fi 



+ 



i - gy P 

F 2 
°p 



2g v 



F 2 



P 2 

-m N 



5 2 
— a„r 

2 



(66) 
(67) 

(68) 
(69) 

(70) 

(71) 
(72) 



To summarize the essential observations: 



• While (gvp — 1) = coincides with the DLFP, (gvui — 1) does not "run" at one-loop 
order. There are two reasons for this. First, there are no contributions to RGE from 
one-loop vertex corrections to the w-NN coupling because all the divergences in wNN 
3-point functions are canceled by those in the nucleon self-energy diagrams figuring 
in wavefunction renormalization shown in Fig. id. Second, there are no uj couplings 
to other mesons at tree order of the ^BHLS Lagrangian, so there cannot be meson- 
loop contributions. However, as mentioned, at higher order this is no longer true. For 

3 This result resembles the one-loop result in QED. In QED with U(l) gauge symmetry, all divergent terms 
of the photon-electron-electron three point function are canceled by the divergent terms in the electron self 
energy diagrams. This cancelation is due to U(l) gauge invariance. However this analogy does not extend to 
higher orders as mentioned in the text. 
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instance there can be uipir one-loop contribution to the w-nucleon coupling involving 
the homogeneous Wess-Zumino (hWZ) term in the anomalous part of the Lagrangian. 

However the hWZ term goes as ~ ( 471 ?p 1 > so the one- loop contribution involving this 
vertex will correspond to normal two-loop order. Therefore we expect {gvu — 1) to run 
slowly, if at all. It may be more appropriate to characterize it as "walking." 

• The RGE for nucleon mass, has the fixed point m^r = 0. But, if T m < at 

the fixed point, = cannot be an infrared fixed point. Note that the w-nucleon 
coupling (gvoj — 1) &j will affect the sign of T m < 0. Since (gvw — 1) does not run at 
one-loop level, T m could become negative at the fixed point unless g w = at that point. 
On the other hand, at two-loop or higher, {gy^ — 1) could be "walking" to zero moving 
toward the DLFP, thereby making JF m > near the fixed point, in which case m,M = 
will become an infrared fixed point. This indicates an intricate interplay between the 
nucleon mass and the w-nuclear dynamics. This is an important point to which we will 
return in Section 15.31 



Figure 1: The diagrams for RGE for u-NN coupling. The thin and thick solid line present the 
nucleon and the background field of the oj meson in the background field gauge respectively. 
The diagrams (a)-(e) contribute to the wave function renormalization of the nucleon field. In 
the calculation of u-NN coupling, the divergences given by the diagrams (f)-(j) are canceled 
by the divergences of the wave function renormalization of the nucleon field. 



4.2 RGEs in the Mirror Assignment 

We now turn to nucleons in the mirror assignment, the parity-doublet model with too 7^ 0. 
Here the treatment is a bit subtler because of non-zero too, which is chirally invariant and 
associated with the explicit breaking of scale invariance tied to the QCD trace anomaly. If 
too <C Atotv where Amyv is the dynamically generated mass, then one could treat too as a 
perturbation. In this case, we will have a qualitatively similar result to that of the baryon in 
the standard assignment. If on the other hand too 3> Atojv, one could adopt heavy-baryon 
chiral perturbation theory. Both cases were discussed in [14] for U(2) symmetry. In between 
the two extreme regimes, we have no reliable tools to address the problem. In this paper, 





+ 



A 
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we are interested in the case where mo is "big'fl since that is what seems to be indicated in 
Nature. 

As in [13] , we will work to the leading order in 1/ thq . With the heavy-baryon field defined 



as 



B + 



exp [irriQV ■ x\ 



N + 



(73) 



the Lagrangian (j37|) takes the form 

C N = iBv^D^B - Am + B + B + - Am_B_5_ 



+ g Vp Bv»a u B + g Vu) Bv» - gj±\ B 

+ g A B Us^ps tanh 6 + v" Pl ol^B , 



where is the spin operator and 



Am± = mjv ± — rriQ . 



(74) 



(75) 



Making slight changes in [13] for the local SU{2) x U(l) symmetry, we can easily write down 
the RGEs: 



dM D 
dp 

dM s 
dp 
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1 

m 



p-r- (gvcu - i) 

dp 







(76) 
(77) 

(78) 

«9) 
(80) 



4 In the vicinity of the chiral restoration point n c , heavy-baryon ChPT could be applied even if mo were 
not so "big." Since we are approaching the DLFP which is supposed to be close to n c , our treatment should 
be reliable. 
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where M 



S/D 



(mjv + ± ?tijv_ ) 2 /4 and and £?i are given by 



^2 
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G-2 



a p9A 



9v P 



F 2 



+ 



' 9\9vp_ 



4 - ^-9v P + 3gv P 



2F3 



P 2 



a 2 + ^ 
9a~\ 
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1 

p2 



+ 25Vp , 

1 - 9V P 
F 2 

crp 



4gy P 
F 2 

crp 



3(l-g Vp ) + 



(81) 



(82) 
(83) 

(84) 



(85) 
(86) 

(87) 

which is the same as what was found in p3]. This follows because of the absence of w-nucleon 
interaction terms for RGEs at one-loop order. From the RGE of {gvu — 1) in (|8U|) follows 
that the coupling (gvui — 1) does not scale at one loop in this case also. Similarly to the 
standard assignment, (gvu — 1) will start scaling only at two-loop and higher-loop order. 



A quick glance at the RGEs in Eqs. (|77p -([80p gives us the infrared fixed point 
{M S ,M D ,1- g Vp ,g A - g Vp ) = (mg,0,0,0) 



5 EoS in the Vicinity of the Dilat on- Limit Fixed Point 

Although, as stated in Introduction, the chiral restoration point n c may be inaccessible 
by local field theory, an interesting theoretic issue is what one can expect at very near 
the DLFP which is presumably close to the VM of hidden local symmetry. This could 
have phenomenological relevance because away from the DLFP, there could be precursor 
phenomena for which the treatment made in this paper is applicable. 

Two issues that can be addressed are the nuclear symmetry energy of asymmetric nuclear 
matter and the repulsion induced by oj exchanges in dense matter. We briefly discuss these 
matters. 

5.1 Symmetry Energy 

The nuclear symmetry E sym figures in the energy per baryon of asymmetric nuclear matter 
as 

E(n, 5) = E (n) + E sym 5 2 + ■■■ (88) 
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where 5 = (N — P)/{N + P) with P{N) standing for the number of protons (neutrons) and 
the ellipsis stands for terms 0{5 ) and higher. We will treat the dBHLS Lagrangian in the 
mean-field (MF) approximation along the line developed in [HI [9]. Since the results are more 
or less the same for both assignments, we will treat only the case of standard assignment. 

As mentioned before and explained in detail in [8j [9] , the MF approximation with chiral 
Lagrangians is justified up to - and slightly above - the nuclear saturation density which can 
be identified as the Landau Fermi-liquid fixed point. Given that we have no reliable guidance, 
we will simply extrapolate it to higher densities and see what can go wrong. From the mean 
field equations of motion for p^, ojq and Xi one can immediately write down the potential 
energy contribution to the symmetry energ}H 



£pot 



9pNN _ (dvp 1) 



n = — —= rrrTl 



^ 2m? '• 2(F* p x*) 2 
where we used 

9* p nn = idvp - l)0p , (90) 
™* = F; p x*g*p (91) 

with x* = (x)V-^x- Here and in what follows, the asterisk denotes density dependence. Note 
that both the numerator and the denominator of ([890 are sliding with density. 

We should point out the difference between the mean field performed here (that we will 
call "xMF" ) and the relativistic mean field found in the literature (that we will call "standard 
MF"). In "realistic" relativistic standard MF approaches, one adds higher dimension field 
operators based on "naturalness conditions" . With the parameters of the Lagrangian taken 
unsealing in density, without higher-dimension fields, the system at normal matter density is 
much too stiff, with a compression modulus more than 3 times what is indicated by nature. 
Here the mean field of higher dimension field operators is to capture the physics of many- 
body forces and interactions and effectively introduces density dependence similar to what 
one has in xMF. See for instance [16]. The higher-dimension field operators in the standard 
MF approaches modify the mass parameter (|91|) with additional density-dependent terms. 
But the shortcoming of these approaches is that there is no constraint that controls the EoS 
at higher densities. 

In our xMF formulation with dBHLS Lagrangian also, we could in principle bring in 
higher-dimension field operators. However here, hidden local symmetry constrains the forms 
of the higher-dimension operators to include, so there is no arbitrariness as in standard MF. 
As recently stressed [1] , the mean-field approach with hidden local symmetry possesses correct 
thermodynamic constraints that allow a controlled extrapolation to higher density, such as, 
for example, the rearrangement terms that reflect many-body terms, including certain many- 
body forces, that are captured in the high-dimension fields in the standard MF approach [9]. 

The distinctive feature of the mean-field result with the dBHLS Lagrangian, Eq. (|89|) . if 
taken seriously at its face value, would imply that as density moves towards the DLFP, the 



5 It is generally agreed that the kinetic energy contribution is small compared with the potential energy 
term, particularly in the presence of strong tensor forces. We will not consider the kinetic energy term in this 
discussion. 
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symmetry energy should strongly decrease, going to zero at the fixed point. Such a softened 
symmetry energy, including the "supersoft" one that drops to zero at ~ 3no, is one of the 
variety of scenarios found with suitable parameters picked in phenomenological models |17j . 
We should mention that although somewhat extreme and perhaps exotic, there is at present 
nothing that rules out such a soft symmetry energy. 

An alternative scenario that is suggested by our results is that the expression (189j) for the 
symmetry energy simply breaks down before reaching the DLFP. When dense matter is de- 
scribed by having skyrmions embedded in a crystal lattice, the skyrmion matter representing 
normal baryonic matter transforms at a density n-yi2 ~ (1-3 — 2)no to a matter consisting 
of half-skyrmions. It involves the transition from a state {qq) ^ to a state with (qq) that 
goes to zero on the average in the unit cell. In the latter, the quark condensate is however 
non-zero locally, so there can be scalar density wave. Chiral symmetry is broken, however, 
with pions propagating, but with a higher-dimension order parameter. 

This skyrmion-half-skyrmion transition at n ly / 2 turns out to have a drastic effect on the 
symmetry energy: The symmetry energy decreases in going toward as does the MF ()89|) . 
but then turns over at and increases strongly, making the E sym stiffer. This signals 
drastic departure from ([89]) . 

It is easy to understand this "cusp" structure in nuclear EFT with the dBHLS Lagrangian 
with the (gvp — 1) coupling going toward the DLFP. In nuclear EFT approach, the nuclear 
symmetry energy is dominated by the nuclear tensor forces that receive major contributions 
from one-pion exchange and one-p exchange. It is roughly given by 

E^ m ~ M (92) 

where Vp is the radial part of the tensor potential and E « 200 MeV is the average excitation 
energy of the states dominantly excited by the tensor forces. The pion contribution and the 
p contribution to the tensor forces, appearing with a different sign, tend to cancel as density 
increases. Thus we can see from (|92|) that the symmetry energy will decrease as the p tensor 
gets enhanced by the dropping p mass. However at n 1/ / 2 , the p tensor coupling gets suppressed 
by the decreas of the (gvp — 1) coupling, thereby leaving the pion tensor fully active above 
n 1 /2- I n @]) this feature was exploited in describing the EoS of compact stars and succeeding 
to explain ~ 2 solar mass neutron stars. The significance of this result is that the mean-field 
result of Eq. (|89p cannot be extended to higher density when there is a topology change of 
the type we encounter, which is not encoded in the standard mean-field approaches. 



5.2 Nuclear Repulsion 

Contrary to the prediction made in [3] with U(2) symmetry for the p and u, the nuclear 
phenomenology of [4] shows that the strong suppression of the w-induced repulsion cannot 
take place as one approaches the DLFP. One can see this also in the xMF with the dBHLS 
Lagrangian. The potential energy per baryon given by an u exchange between two nucleons 
is of the form 

E p 0t = |w! n + . . . (93) 
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where the ellipsis stands for three-body and higher-body contributions which will be of 

* 2 

( 9 $ N »i ) or higher order. Using the tree-order mass formula for the in-medium uj, m^ 2 = 
(gl F Lx*? and the coupling g* NN = gZidy^ - 1), we have that 



+ (94) 



We expect that while the numerator will drop only slowly - if at all, the denominator will 
drop somewhat less slowly. Thus one sees that the u repulsion will not undergo an appreciable 
change as the DLFP is approached. 



5.3 Interplay between the Nucleon Mass and the Repulsive Core 

Although we have neither mathematically rigorous arguments nor any results of phenomeno- 
logical studies using the formulation of [4], we can still make a qualitative connection between 
the main source (or origin of ) of nucleon mass and the property of the w-nucleon coupling in 
medium. As discussed recently [18j . there are at least two reasons to think that the nucleon 
mass has a large component that is chiral-scalar and hence does not disappear at the chiral 
transition. One is a lattice result: By artificially unbreaking chiral symmetry in a dynamical 
lattice simulation [19] (see also [20]), Glozman et al. find that baryons - and also mesons - 
remain massive after chiral symmetry is supposedly restored [21]. In fact, phrased in the form 
of Eq. (I26p . itlq is found to be surprisingly large. The other is evidence from dense skyrmion 
matter where one notes that as density increases, the energy per baryon of the system scales 
with the effective pion decay constant that falls quite slowly once the system enters into the 
half-skrymion phase [18] . 

As proposed in [U [18] and stated above, the presence of a large chirally invariant term 
in the nucleon mass can be effectively captured in the parity-doublet baryon structure with 
the nucleon mass given by f)26|) . tun = tuq + Am(rt). In principle, one could determine mo 
from low-energy hadronic processes involving parity-doublet baryons with (IBHLS. For this, 
high-order chiral perturbation theory with the vector mesons suitably incorporated would be 
needed. Such a formulation is not yet available. 

A rough range available up to date is that tuq could range from ~ 200 MeV to ~ 800 
MeV [31CII1 [22]. 

Let us briefly discuss what we can say from nuclear phenomenology with, say, tuq ~ 800 
MeV. In [1] , the structure of dense skyrmion matter simulated on crystal lattice treated with 
a truncated form of dBHLS [2] was translated into the scaling behavior of the parameters 
figuring in the dBHLS Lagranigian. The rationale for such a procedure was that the topology 
change from skyrmions to half-skyrmions present in the crystal matter could be interpreted 
as a change in the parameters as density increases across ni/ 2 - For the nucleon mass, the 
scaling was taken to follow the empirical trend up to the normal matter density no, giving 
m* N /m,N ~ 0.8 at n = no, and then assumed to remain unchanged after n\/2 ~ 1.2uq. We can 
consider this choice as corresponding to mo/rriN ~ 0.8 in the parity-doublet dBHLS model 
with Am — > as n — > n\ii at which (qq)* vanishes^. For this value of m* N - and with other 

6 Although the condensate {qq) is not a proper order parameter for chiral symmetry for the crystal matter, 
we will identify it to be associated with the dynamically generated mass term Am 
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parameters fixed, it was found that even a slight decrease in the w-nucleon coupling would 
generate instability of nuclear matter for n > n\i2- This is an indication that for the given 
value of m* N , the g^NN cannot run much in density after n\/2- How this interplay between 
the nucleon mass and the w-nuclear coupling manifests in nuclear dynamics can be seen in 
RGEs and is described schematically in Appendix. 

The reason for this phenomenon is easy to understand. In nonrelativistic description in 
nuclear EFT, lowering m* N engenders strong suppression in attraction, the suppression factor 
going like ~ (m* N /tun) 2 ■, basically a relativistic effect on the scalar density at higher baryon 
density. On the other hand, lowering g^ NN decreases, roughly quadratically, the w-exchange 
repulsion, causing instability of the matter at high density. Therefore with other parameters 
fixed, the two are tightly correlated in the EoS. Given that mo is not known at present, 
one may vary both mo and g^NN 111 sucn a wa Y that the EoS remains within the ranges 
provided by heavy-ion experimental data available up to ~ 4no and consistent with neutron- 
star observables. This could offer some indication as to whether {gvu — 1) runs or rather 
walks as we found in this paper. This feature is illustrated by a simple mean-field treatment 
with dBHLS in Appendix. It should be emphasized that the mo is the part of mass that has 
nothing direct to do with the spontaneous breaking of chiral symmetry, so if it is as big as 
80% of the nucleon mass, the origin of the bulk of proton mass is yet to be explained [18] . 
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Appendix 



In this appendix, we examine in more detail, with the help of the RGE Eq. (|62p . the possible 
interplay between the nucleon mass and the w-nucleon coupling. For simplicity we consider 
the standard scenario with mo = 0, but it is straightforward to extend it to the case with 
mo 7^ 0. In doing this we will continue to assume that DLFP is an infrared fixed point. It 
cannot be an ultraviolet fixed point because 
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(95) 
(96) 



does not become negative at DLFP unless (gvu — 1) 2 <?S goes to oo at DLFP. 
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Figure 2: The nucleon mass vs. w-nucleon coupling as a function of density for varying 
density dependence of gy u - What is notable is that the nucleon mass stops dropping at a 
density slightly above nuclear matter density uq and stays more or less constant above that 
density. The density n^A) corresponds to h referred to in the text. 



Let us first consider the case where the cj-nucleon coupling is non-scaling , i.e., a constant 
or B = in the parametrizatioijl] defined in Fig. [21 

9 V" ~ 1 = 1 (Q7) 

g° Vu} -l l + Bn/n { } 

where the asterisk stands as before for density dependence. As the infrared fixed point, 
rriM = 0, is approached from low density, T m will remain positive until one reaches a certain 
point, say, 'A' with density n(A)- The density ni^) cannot be too low because \j? which should 
correspond to the momentum cutoff A 2 should be big compared with the other terms with 
negative sign. Up to nrj^i m N will decrease. However, beyond A, JF m will become negative, 
and tun = will no longer be an infrared fixed point since 

dm N 3mjv ^ fnQ , 
^=16^ <0 - (98) 

This means that from the point A, the nucleon mass either increases^ or more or less stays 
constant because the sign of fJ-^rf- after passing the point A is opposite to the sign of fJ-^jf^ 
before passing the point A. 

7 This is the possibility considered in the nuclear effective field theory calculation of EoS for compact-star 
matter in [J with ri(A) identified as the density at which the half-skyrmion phase sets in in dense skyrmion 
matter. Note that this parametrization would not be valid for large densities. 

8 This increase may be an artifact of the approximations and most likely unphysical. This point is being 
further studied. 
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Now consider the case where the w-nucleon coupling decreases. Up to n^, the nucleon 
mass will decrease as before. But, passed the point A, there are two possibilities depending on 
the speed with which the w-nucleon coupling decreases. If the w-nucleon coupling decreases 
'fast' enouglJlto have - and keep - J- m positive, then the nucleon mass will continue decreasing 
after the point A. But, it will decrease slowly because T m will be small. If, on the contrary, 
the w-nucleon coupling decreases very slowly, which will make T m negative at some point C 
with U(n\ > ntM, then the nucleon mass will not decrease after the point A. 

Because of the interplay between the nucleon mass and the w-nucleon coupling given by 
RGEs, the slowly decreasing w-nucleon coupling makes the nucleon mass also decrease slowly. 
This causes the DLFP to be delayed and reached at higher density, because the decrease of 
gv P and that of gA, given, respectively, by (f63|) and (fM|) . are related to the decrease of the 
nucleon mass, Therefore, the p- nucleon coupling will decrease fast as in [31 E] before reaching 
n( A ), the density corresponding to the point A, but after passing n^, it will start to decrease 
slowly. It follows that the DLFP gy p — 1 = with broken U{2) symmetry will be reached at 
higher density than the density for DLFP with U(2) symmetry. 

This subtle interplay between the nucleon mass and the w-nucleon coupling can also 
be seen in the mean-field calculation with dBHLS Lagrangian of Section [3] in the mirror 
assignment with itlq. The nucleon mass is connected to the w-nucleon coupling by the equation 
of the motion for x an d w, and the in- medium property of the x condensate - equivalently the 
in-medium mass of the dilaton - controls the behavior of the nucleon mass at high density. 
The nucleon mass depends on x — (x) as given in (|58l) : 



m N± ==F 5 r 2 p L x + y (gi^-xj + m l- (") 

In the mean field approach, taking the limit tr [SS T ] — > is replaced by x ~ > 0- The slowly 
decreasing w-nucleon coupling not only causes the nucleon mass to drop slowly, but also 
delays the dilaton limit, gA = gvp, to higher density. 

In Fig. [2] is shown a simple illustration of the above observation. For a given B, the 
nucleon mass is calculated by fitting the binding energy and the pressure of nuclear matter 
at no0 It is remarkable that the nucleon mass drops almost linearly in density up to slightly 
above n ~ no, say, n^j^j, and then stays more or less constant. In terms of the parity-doublet 
nucleon mass defined by Eq. (|26|) . we have mo/wiN ~ 0.7. The interplay between the nucleon 
mass and the w-nucleon coupling as revealed in this way is uncannily similar to what was found 
by the RGE analysis given in this paper and consistent with what was phenomenologically 
observed in nuclear EFT description with a BR scaling modified by topology change [1]. 

9 Here, 'fast' doe not mean that it decreases faster than the p-nucleon coupling which is to drop quite fast. 

10 The two density-dependent quantities involved are m* and gy w that are determined by the binding energy 
and the pressure at n = no for given B. The behavior of the nucleon mass of Fig. [2]then follows. 

n In the text, this was referred to as h. We conjectured there that it corresponds to the skyrmion-half- 
skyrmion transition density rii/2 
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